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Abstract 

We derive new crystal melting models from Chern-Simons theory on the three-sphere. 
Via large N duality, these models compute amplitudes for A-model on the resolved 
conifold. The crystal is bounded by two walls whose distance corresponds to the 
Kahler modulus of the geometry. An interesting phenomenon is found where the Kahler 
modulus is shifted by the presence of non-compact D-branes. We also discuss the idea 
of using the crystal models as means of proving more general large N dualities to all 
order in g s . 



1 Introduction 



Topological string theory [H [2] is currently undergoing a drastic paradigm change. 
Reshetikhin, Okounkov, and Vafa [5] realized that various amplitudes for topological 
A model on C 3 can be expressed in terms of classical statistical models of a melting 
crystal. Iqbal, Nekrasov, and Vafa |4] proposed to interpret the crystals in terms of 
quantum foam or Kahler gravity, which is the target space theory of A-model closed 
string theory. Mathematically speaking, this means that Gromov-Witten invariants 
are related to the so-called Donaldson-Thomas invariants j5j E] . 

Central to the dramatic paradigm shift in topological string is the interpretation 
of the Calabi-Yau crystal as describing the violent fluctuations of topology and the 
geometry at microscopic scales. This is reminiscent of geometric transition, where 
open string theory and closed string are related via a local change of topology and 
geometry. Or rather, when crystal picture is combined with geometric transition, one 
naturally expects that the geometric change is part of the gravitational fluctuations or 
quantum foam. In this paper, we realize this expectation and make it precise. 

We propose a crystal melting model that describes the A model closed strings on 
the resolved conifold O(-l) © 0( — 1) — > P 1 . Our model is a simple modification of 
the model for C 3 . The Kahler gravity interpretation leads one to view the gluing pre- 
scription of the topological vertices as computing partition functions of crystal models 
for general toric Calabi-Yau manifolds. Although the resolved conifold was discussed 
in that context in [3], what we propose is different from the prescription described 
there. Indeed it was emphasized in [7] that "the global rule of melting is absent for 
closed strings on toric Calabi-Yau manifolds with more than one fixed point of the toric 
action" . One of the purposes of the present paper is to amend this situation. 

Our model is obtained from the large N dual Chern-Simons theory on S 3 . We show 
that the Chern-Simons theory can be formulated as a simple unitary matrix model that 
involves a theta function. This representation is then used to obtain a free field formula 
for the Chern-Simons theory, which is interpreted in terms of a statistical model. 

It is also possible to introduce non-compact D-branes to the crystal, enlarging the 
arena of study to include open strings. It was shown in [7J for C 3 that this corresponds 
to having defects in the crystal. In Chern-Simons theory, the observables are the Wilson 
loops that go around the circles in various knots and links in the three-manifold. We 
show that the computation of a Wilson loop along an unknot can be nicely done in the 
unitary matrix model. This then translates to a natural crystal model with defects that 
represents some number of non-compact D-branes intersecting the P 1 in the resolved 
conifold. The fact that these D-branes fit neatly into the crystal shows that our model 
of the Calabi-Yau crystal is a natural one. 

The crystal melting model is also a useful computational tool. For one crystal model, 
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there are two ways to represent it in terms of free bosons and fermions. We use this 
freedom to explicitly compute certain amplitudes in Chern-Simons theory. We find an 
interesting phenomenon where the Kahler modulus of the resolved conifold is shifted by 
a multiple of g s in the presence of non-compact D-branes. We also discuss the possible 
application of the crystal representation to prove more general examples of topological 
string large N duality to all order in g s . 

The plan of the paper is as follows. In section [2] we propose a crystal melting model 
and demonstrate that it computes the partition function for the resolved conifold. In 
section [3] we explain how the crystal picture naturally arises from the dual open string 
theory. We also discuss the non-perturbative mismatch. In section 0] we derive from 
Chern-Simons theory the crystal models for non-compact D-branes, realizing them 
as defects in the crystal. Section [5] discusses the possible application of the crystal 
computation as a way to prove large N dualities to all order in g s . 

2 Crystal melting model for the resolved conifold 

Let us recall the crystal model for C 3 [3] . The zero-energy configuration is the positive 
octant x, y, z > in R 3 filled with atoms. Here an atom at (xq, yo, zq) is a filled box 
{(xo + s x , yo + s y , z + s z )\0 < s x ,s y ,s z < 1}. We consider removing atoms from the 
corner. The allowed configurations are defined recursively as follows: The configuration 
where the whole octant is filled is allowed. If an allowed configuration has an atom at 
(xo, yo, Zo) such that there are no atoms in the region {(x, y, z)\x < xq, y < yo, z < zo], 
one can remove the atom at (xo,yo,Zo) to obtain another allowed configuration. The 
allowed configurations are also called 3D Young diagrams, in analogy with the familiar 
counterpart in two dimensions. The partition function is 



where the summation is over 3D Young diagrams ir, and q = e~ 9 % \ir\ is the number of 
atoms removed. This partition function agrees with the partition function of A-model 
closed strings on C 3 . This fact can be proved by the use of free field techniques familiar 
in string theory [3]. Below we generalize the technique to the situations of our interest. 

The model we propose for the resolved conifold is the following. We add one more 
condition that further restricts the allowed configurations: Atoms in the region x > N 
cannot be removed. Here N is related to the Kahler modulus t as t = g s N. Note that 
this condition introduces a "wall" that together with the original three walls constitutes 
the toric diagram for the resolved conifold. 

Now we demonstrate that this crystal melting program indeed reproduces the par- 
tition function for the resolved conifold. For this purpose, we express the partition 




(2.1) 



7T 



2 





(b) 

Figure 1: (a) The crystal melting model for the resolved conifold 0(— 1)©0(— 1) — > P 1 . 
The edges, drawn as solid lines, of the positive octant bounded by the wall at y = N 
form the toric diagram of the resolved conifold. (b) Many atoms have been removed 
from the crystal. Atoms cannot be removed from the region beyond the wall at y = N. 

function in terms of free fermions and bosons: 



reZ+1/2 r 
{?/V,^ s } = S r+Sj0 , 

<p(z) = x — ia log z + i 



«n 

nz r 



[a m , a n ] = m5 m+nfi . 



(2.2) 

(2.3) 
(2.4) 

(2.5) 



These are related via 



id(p(z) =: ^{z) :, tf}{z) =: e irp(z) :, ijj(z) =: e" 



■icf>(z) 



(2.6) 



Now we define 



T±(z) = exp y a± n . 



n>0 



(2.7) 



It is well known that neutral (zero momentum in the bosonic language) fermionic Fock 
states are labelled by (2D) Young diagrams \i which we denote as /i = (//i > //2--- > 
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fid > 0). More explicitly, such Fock states are given by 

oo 

i/i) = n^-w-i/*i°» 



'=i 

d 



i=l 



where |0)) is the state that is annihilated by all ip r , r G Z + 1/2, and we have defined 

Oi = iH-i + l/2, hi = n\ - % + 1/2. (2.9) 

/x* is the transposed Young diagram. The Virasoro zero mode Lq counts the number 
\li | of boxes in the Young diagram \x: 

A)|m) = I*}- (2.io) 

Two Young diagrams A and /j are said to interlace (and we write A >- /i) if they 
satisfy 

Ai > /xi > A 2 > /i 2 > ... (2.11) 

In other words, A and /i interlace if and only if A contains \i and /i contains A with 
the first row removed. The interlacing condition is equivalent to the local condition 
for two Young diagrams one finds by slicing the allowed configuration of the crystal by 
the planes x = y + j and x = y + j + 1 [3]. The operators T±(z) are useful because of 
the properties 



r + (i)|A> = $>>, 

r_(l)|A> = 5>). (2.12) 
The partition function for the crystal can be written as 

Z crystal (q,t = g s N) = (0| fn 9 L »r + (l)] q L ° ( T_(l)q L ° ] |0) 



v,n=l / \m=l 



N 



(oin r +(^ i/2 ) n r ^~ (m ~ i/2) )i°>- ( 2 - 13 ) 



n=l m=l 



This can be understood as slicing the crystal by planes x = y + j, j G Z. Note that 
we have a finite product of vertex operators acting on |0). This restricts a 3D Young 
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Figure 2: (a) Closed string slicing: This slicing by planes x = y + j allows one to 
compute the closed string amplitude as in eq. (I2.14p . (b)Open string slicing: Another 
slicing of the crystal by planes x = z + j corresponds to the free field representation 
eq. (13.1 5p obtained from Chern-Simons theory. 

diagram to have a trivial 2D Young diagram on the slice x = y — N. The interlacing 
conditions then imply that the 2D Young diagrams must have at most one row on the 
slice x = y — N + 1, two rows on x = y — N + 2, etc. Thus the free field correlator 
represents a crystal model bounded by a wall at y = N. See figure [2](a). 
Now we can explicitly compute the partition function. 

a n V"^ ± — Q 

Z CIysta x(q,t = 9sN) = (0|e~^>°^Me-^> o ^r Q -"|0) 

r 1 -" Nn 

= Q^ n>0 n[n] 2 

= M(q)e ^">°^F. (2.14) 

Taking iV — > oo pushes the wall at y = N to infinity, and the partition function reduces 
to the result for C 3 . 

This crystal model and the resulting amplitude are different from those discussed in 
[4j. While our crystal has a fixed finite size in the y direction, in [4] the distance between 
the two crystal corners are not fixed because two finite size 3D partitions are connected 
through a region of length t = g s N. Consequently, in stead of a single power of M(q) in 
our model, the model [4J gives the square of M(q). More generally, a closed string par- 
tition function contains M(q) x ^ x ^ 2 , where xPO is the Euler characteristic of the target 
space X. If the target space X is non-compact, the definition of the Euler characteristic 
is ambiguous. In the context of large iV duality, it is known [8] that one should assign 
the value 2 to the Euler characteristic of the resolved conifold as we just did. This is nat- 
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ural in the sense that the target space admits one Kahler deformation but no complex 
structure deformation, and the general formula for a (compact) Calabi-Yau manifold 
is xPO = 2 [(#Kahler deformations) — (^Complex structure deformations)]. 

3 Large N dual open string theory 

In this section, we study the crystal melting problem from the point of view of large 
N duality 

The crystal model in the previous section can also be expressed as 

^crystal (q,t = g s N) 

oo oo 

= (0|nr + (g n - 1/2 )i^< J vn r -(^ (m_1/2) )l°) ( 3 - 15 ) 

n=l m=l 

Here I^kn is the operator that projects onto the subspace spanned by \fi) such that 
the Young diagram \i has at most N columns. This free field expression corresponds 
to slicing the crystal by planes z — x = j, j G Z. See figure 12(b) . 

We call this the "open string slicing" because, as we will see below, this representation 
of the crystal naturally arises from Chern-Simons theory. 

3.1 Unitary matrix model for Chern-Simons theory 

The large N duality of Gopakumar and Vafa relates U(N) Chern-Simons theory on 
S 3 to topological closed string on the resolved conifold 0(—l) © 0(—l) — > P 1 . The 
dictionary is that the Kahler modulus t of the closed string theory geometry is identified 
with the 't hooft parameter g s N. Certain amplitudes on the resolved conifold, including 
the closed string amplitudes, can be computed within the framework of Chern-Simons 
theory. We now develop a unitary matrix model formulation of Chern-Simons theory, 
which will be used to derive the crystal model for the resolved conifold later. 
The partition function of the U(N) Chern-Simons theory on S 3 is given by [9] 

Z cs (N,k,U(N)) = {k + x )N/2 U^^- (3-16) 

Here 2iri/(k + N) = g s is the string coupling constant, and the product is over the 
positive roots of SU(N) C U(N) which are given by o;^ = — tj G C iV (~ Cartan 
subalgebra) for i < j. p = (1/2) J2 a >o a = J2iLi(.^~2^ ~ ^) e * ^ s ^ e Weyl vector. Now 
note the following formula by Weyl for the denominator of the Lie algebra characters 

Y[ 2 sinh(a -u)=J2 <™)e w{p) ' u , (3.17) 

a>0 weW 
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where W is the Weyl group isomorphic to the permutation group Sn and u is an 
arbitrary element of the Cartan subalgebra of U(N). We use this identity to rewrite 
Z cs {N,k,U{N))=:Z cs {g s ,N) as 



iV(jV-l)7ri 

e " 4 



e[w ) e 9Mp>P 



Zcs(9s,N) - J2 

= (g) e~%--5-^ cs Gfc,iV) (3.18) 
Here we have factored out the non-trivial part of the partition function: 

Z cs (g s ,N) = e(w)q^ w ^ 2 . (3.19) 

q is again e~ 9a . In what follows, we "analytically continue" in g s and regard q s as a 
complex parameter with a positive real part. We can introduce another sum over the 
Weyl group and an integral over the maximal torus as follows^] 

Zcs(g.,N) 



W 



w,w'£W 
N 

I , I IS \ 

J(w(p)-w'(p))-< 



1 1 J i=l V 7 w.w'eW 



(3.20) 



Here 



,2 



^oo(e^;g):=^g^ me (3.21) 

meZ 

is one of Jacobi's theta functions. 

By making use of the Weyl denominator formula eq. (13.171) again, we get 

N 



*» = ^/(n§^*))(n2-^) 2 - (3.22) 



w\ 



vi=l / \a>0 



The second factor now represents the Haar measure for U(N) pushed down to the 
maximal torus. The partition function can be written in a very simple form 



^cs 



/ dUdet$ 00 {U;q), (3.23) 

JU{N) 



1 Here we use the identity q m I 2 — J^ n |f ^oo(e* e ; q)e" nS . If we instead use 
q m I 2 — J^ g e~2a7e m ", we get the matrix model with a non-compact integration region intro- 
duced in [10[ 111] . The matrix model there can related be transformed to our unitary matrix model 
via u = i{9 + 27m), performing the sum over nGZ and a modular transformation. The author thanks 
Hirosi Ooguri for pointing this out. 
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where the measure is normalized so that the volume of U(N) is unity. This expression 
holds for any gauge group of the Chern-Simons theory on S 3 when the corresponding 
Haar measure is used. 



3.2 Crystal from Chern-Simons theory 

Now we use the product formula for the theta function 

oo 

tfoo(e ie ;g) = n( 1 -9 i )( 1 + ei V'" 1/2 )(i + e- i V~ 1 



JI(l-^)]exp 



D- 1 : 



^in6 _|_ g— ind 



n>0 



n\n\ 



(3.24) 



to write 



^cs 



VJ = 1 

' oo 

m 

K j=i 



-q J 



N 



A' 



dU exp 



dU exp 



,n>0 



ran 



£ 

.n>0 



Trf/ n + Tt jj- 



n\n\ 



(3.25) 



To obtain the free-field expression for the partition function, we introduce the co- 
herent states: 



\U) := exp 



V -Tr U n a. 
^ n 

n>0 



|0). 



These states satisfy 



(3.26) 

(3.27) 
(3.28) 

where ld<N is the projection to the subspace spanned by such that the number of 
rows in /i is less than or equal to N. This formalism was extensively used in the context 
of 2D Yang-Mills theory which has recently been attracting some attention p^2l [13l fl^] . 
See [15] and the references therein. 
By making use of \U), we can write 



a n \U)=TrU n \U), 
dU\U)(U\ = W, 



N 



Y[{^--q j )) (0|e E ">°^T / c/[/|[/)(f/|e^" >o ^r|0) 



1 £- — ' Ti 



N 



- q j ) ) (0|e E ">° »W l d < N e^ n>0 »M |0) 



(3.29) 
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Since a n — > —a n is equivalent to R — > R* (c.f. eq. (12.81) ). we finally obtain 

/ 00 \ N 

Z cs = \Y[(l-<l j )) (0|e~ E ">°5il dt < iV e" E "> o ^r|0) 

= Z(q)- N Z CTystal ( q] t = g s N). (3.30) 

We have demonstrated that the open string slicing eq. (13. 15ft naturally arises from 
Chern-Simons theory. Note that there is a mismatch by the factor £,(q)~ N between 
Zqs and Crystal, where £(g) = l/n^Li(l — * s the "renormalization factor" which 
was found in [7j to be associated with a non-compact D-brane. As in [7j, we use the 
modular property of 77(g) = g 1 ' 24 ^?) -1 , namely 77(g) = y^2n/ g s r](q), q = e~ 47r l g % to 
argue that it does not contribute to the perturbative amplitudes at genus no less than 2 
when comparing the open and closed string sides. For genus amplitudes, the mismatch 
is absorbed into the usual ambiguities. 



4 Adding D-branes 

We can add non-compact D-branes to the system. In the language of Chern-Simons 
gauge theory, this corresponds to placing Wilson lines going through circles of links. 
In the case of an unknot, we will be able to see the connection to the description in [TJ. 

On the open string side, we consider placing a stack of M non-compact D-branes in 
T*S 3 intersecting the S 3 along an unknot S 1 [16J. Since the new D-branes are non- 
compact, we treat them as non-dynamical, acting as a source to the gauge fields on S 3 
via an interaction. This interaction is obtained by integrating out the degrees freedom 
coming from the open strings stretching between the compact D-branes wrapping the 
S 3 and the non-compact D-branes. Let U & U (N) and V £ U (M) be the holonomies 
along the unknot for the gauge fields on the compact and the non-compact D-branes, 
respectively. Then the interaction can be represented as 

The expectation value can be expanded, with the help of Frobenius' formula, as 

<e E ^ Tr ™"> = ^(Tr^)Tr^. (4.32) 

Here Tr M denotes the trace in the representation of U(N) or U(M) specified by the 
Young diagram \x. 

It is natural to expect that (Tr^U) in eq. (I4.32p is computed by the unitary matrix 
model in subsection 13.11 by inserting Tr^U. We now show that this is indeed correct, 
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however with the subtlety that the Wilson line and hence the non-compact D-branes 
have non-canonical framing. 

The object we would like to compute is 



dU det VooiU; q)Tr^U. 



(4.33) 



Going back to the eigenvalue integral, this is 



N 



d9 



n^oo(e^)det [(e^)^] det [(e" 



■iO;\N-i 



., det [( e ^)w+^-'] 



det [(e 



(4.34) 



where we have used the Jacobi- Trudy formula Tr^diag(xi, xn) = s^(xi, xn) = 
detxf+^/det x 1 ^ 1 for the Schur polynomial After cancelling factors between the 
numerator and the denominator, and performing the integrals the matrix integral re- 
duces to 

N 



3=1 



N 



sgn ( xJ]g5(^-i+-0')) 2 

creS N j=l 

det \qh0n-i+3) a 



Up to /^-independent factors, this equals 
g iE£i«(w-«+w) det 



(4.35) 



(4.36) 



The power of q can be written as g^M+^HV 2 , where = 2 i) = S« failM ~ 

2i + 1). This is the factor one obtains when the framing of the Wilson loop is shifted 
by one unit [5J- The determinant is of the form that appears in the numerator of the 
Jacobi-Trudy formula. Hence we have shown that 



/ dU det $ 00 (U; g^i^U 

J dU det &oo{U; q) 
g ( ^ +Ar| ^ l)/2 Tr M diag(g-^ i ,g 



iV-l N-3 
2 



N-1 . 

,q 2 



(4.37) 



Relative to the result for the canonically framed unknot [16J, we see that the matrix 
model computes amplitudes in the framing shifted by one unit. 

This vacuum expectation value of the Wilson loop can be represented as a crystal 
melting model as follows. 



2 A good reference on symmetric functions and the group theory relevant to us is [17]. 
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dUdest $ o(U; qyVtyU 
^(g)^(0|e E ">°S J rft/Tr^f/|f/)(t/|e E ">°^r|0) 

r 1 00 

e(g)^(0|e E "> o S / dt/^— x M (C(fc)) JJa^|t/)(f/|e E ">°^M|0), (4.38) 



where is the S\^\ character of the representation specified by /i, k = (ki, k 2 , ...) is an 
infinite vector with non-negative integer components, and C(k) is the conjugacy class 
of S^i specified by k. Now the powers of a-,- can be moved to the left to act on (0|. 
This yields 

J dUdet# 00 (U; g)Tr M C7 
= ^^"^(/ile^ ^! / c#/|[/)([/|e E ™> o ^M|0) 



£(g) (/i|e^">° ^l d <7ve^™ >0 ^ |0) 

at 4 a n \ ^ a — n 

f (?) <A* l e ^ !rf'<ive~ ^' l>0 |0) 
00 00 

e^r^vi n r +^ n " 1/2 )^<^ n r-r (m - 1/2) )io) 

n=l m=l 



crystal 

where we have defined 

00 00 

z^ir ■= q-^-Mtfl nr + (g n )i^<^ n r-Gr^lo) (4.40) 

n=l m=l 

This free field correlator together with the power of q represents, in the open string 
slicing, the partition function of the crystal melting model whose initial configuration 
is shown in figure [31 The power of q ensures that the initial configuration has zero 
energy. 

It is possible to express the multi-D-brane crystal in the closed string slicing, which 
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Figure 3: The initial configuration of the crystal with defects representing multiple 
non-compact D-branes intersecting P 1 in the resolved conifold. The defects introduce 
faces at y = fi\ = Ni - M + 1, fi f 2 = N 2 - M + 2, ...,fi? M -i = N M -i - 1, a4/ = n m- 

is a slight generalization of the free field representation in [7]. 

yT>— branes 
crystal 

oo Nm 

= <o| \{ r + (g"- 1/2 ) n r_(g-( m - 1 /2))r + (g-(^ + i/2 )) 

n=l m=l 

Nm-i Nm-2 

x ] J r _ (g -(— i/2) )r+(g -(/v M _ 1+ i /2)) |j r_(g-^- 1 / 2 ))r + ( g -^-^ 1 / 2 ))... 

m=N M +2 m=N M -i+2 
Ni N+M 

x Yl r_( g -( w - 1 / 2 ))r + ( g -( JVl+1 /2)) J] ^(g-^- 1 / 2 ))^}. (4.41) 

m=7V 2 +2 m=7Vi+2 

Here /4 = JVi - M + l,/4 = iV 2 - M + 2, ...,fi t M _ 1 = N M -i ~ l,/4f = n m- In the 
closed string slicing, it is possible to explicitly evaluate the correlator to write it as a 
product. This also provides us with an interpretation of Ni as positions of D-branes 
and exhibits an interesting shift in the Kahler modulus: 
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r?D— branes 
crystal 



n n T^^iln n 

«9) M 



~y Qti+m— 1 I J. J. J. J. J Q'm—Ni — 1 

n=l l<m<N+M,m^Nj+l y / i=l iV i +2<m<AT+iW',m^iV,+l 



n ^-^ 

.i<j<i<M 

x «w «p (- e ^ J n «p (e — ^ — J • <" 2 > 

Here we have defined a { := g s (Ni + 1/2), i = 1, M and t := g s (N + M) = t + g s M. 
Again, £(g) can be essentially ignored in the perturbative computation due to the 
modular property of 77(g) = q 1 ^ 2A ^(q)~ 1 - The second factor rii<j(l — e aj " ai ) is also 
present in the multi-brane case of [7j, and written in this way is independent of g s . 
This is the amplitude for M non-compact D-branes in the resolved conifold, which can 
be defined as the Kahler quotient 

{(Xj) e C 4 : |Xi| 2 + \X 2 \ 2 - \X 3 \ 2 - |X 4 | 2 = Re t}/U{l) (4.43) 

with t/(l) action by charges (1, 1, —1, —1). The geometry of the D-branes is [18] 

|X 1 | 2 -Re(a i ) = |X 2 | 2 -Re (i - en) = \X 3 \ 2 = \X 4 \ 2 , ^argX/ = 0. (4.44) 

1 

One thing that is interesting in our computation is that the Kahler parameter is shifted 
from t = g s N to t = t + g s M. It has been known (see, for example, [19]) that the 
presence of D-branes can shift the effective size of the geometry by the string coupling 
times the number of D-branes. Here we have found another such phenomenon. The 
genus zero part of eq. (14.421) in the case of a single D-brane agrees with the results in 

m- 

The fact that that non-compact D-branes can be nicely incorporated to the crystal 
confirms that our crystal model of the resolved conifold is a natural one. 



5 More general large N dualities, instanton count- 
ing, and geometric engineering 

So far we have been discussing the Calabi-Yau crystal in the context of Gopakumar- 
Vafa duality (T*S 3 <^> 0(-l) ®0(-l) -> P 1 ), the simplest example of large N duality 
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in topological string theory. There is a family of generalizations of the large N duality 
which is worth considering in relation to Calabi-Yau crystal. The example of Gopaku- 
mar and Vafa is simple enough to prove the duality (at least at the level of free energies 
and some open string amplitudes) by direct calculations. However, our derivation of 
the resolved conifold crystal from Chern-Simons theory can be viewed as a complicated 
way of proving the duality. In this section we discuss the possible application of the 
ideas in the present paper to prove more general large N dualities. 

Aganagic, Klemm, Marino, and Vafa made a conjecture in [TT] that that the duality 
of Gopakumar and Vafa still holds after taking a Z n orbifold on the both sides of duality. 
On the closed string side, this produces A-type topological closed string theory living 
on the particular fibration of the A n _i ALE space over P 1 . The geometry has n Kahler 
moduli, the sizes of the base P 1 and n additional P 1 that blow up the A n _i singularity. 
On the open string side, we again get Chern-Simons theory, this time living on the lens 
space L(n, 1) ~ S 3 /Z n . Also after taking the orbifold, the relevant open string theory 
is a sector of Chern-Simons theory which contains one classical solution. A classical so- 



lution can be specified by a holonomy exp[27ri/iVdiag(l, 1, 2, 2, n, n)] along 
the generator of the homotopy group. The Kahler parameters are then to be identified 
with linear combinations of the t' Hooft parameters g s Ni,i = 1, ...,n. 

The n = 2 duality was tested via perturbative computations by the authors who 
proposed the duality [TT]. For general n and a related duality, checks have been done 
by showing that the matrix models describing the sector of Chern-Simons theory leads 
to the spectral curves which are the non-trivial parts of the Calabi-Yau manifolds mirror 
to the A- model closed string geometries [211 [221 E3] . The worldsheet derivation of the 
Gopakumar- Vafa duality [21] has also been generalized for these large N dualities |25j . 



There are n + 1 choices (m = 0, 1, ...,n + 1 in the notation of [26J) one can make 

when one fibers the A n _i ALE space over P 1 . The closed string geometry that is 
dual to the S 3 /Z n Chern-Simons theory is precisely the fibration m = [22J that 
was shown to correspond to Nekrasov's instanton counting [23 I2H] for the 5D SU(n) 
gauge theory with vanishing Chern-Simons term [23 H]. (For the correspondence with 
non-zero Chern-Simons term, see [30].) Nekrasov's correspondence between topological 
closed strings and 5D gauge theory has been discussed in [311 [321 S3] by making use 
of the topological vertex [19]. As discussed in the introduction, the computation via 
the topological vertex is closely related to the Calabi-Yau crystal. In particular, the 
computation takes the form of an expansion in q = e~ 9s . 

As we saw in the previous section, the Chern-Simons theory also naturally leads to 
an expansion in q. Hence, it is plausible that one will be able to prove the generalized 
large N dualities to all order in g s by proving that the partition functions are the same 
on the both sides as functions of q [31] . 
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Let us also make an observation on the appearance of the unitary matrix model. In 
[35] , the question of finding matrix models that compute the Seiberg-Witten solutions 
of M = 2 gauge theories was addressed. The matrix models in [11] can be regarded as 
computing amplitudes in the 5D gauge theories with the same number of supercharges. 
By taking a double scaling limit, which is the familiar field theory limit of geometric 
engineering [36], one can compute amplitudes for 4D M = 2 gauge theories from these 
matrix models. By using the technique in this paper, it is possible to rewrite the matrix 
models in [XT] as unitary matrix models. These are similar to, and can be regarded as 
generalizations of, the unitary matrix model (Gross- Witten one plaquette model [37] ) 
that was considered in [35] for the SU (2) gauge theory. 
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